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Some thermodynamic quantities of nonrelativistic ideal boson and fermion gases in the static 
^ , Taub universe are derived to first order in a small anisotropy parameter d which measur- 

ing the deformation from the spherical Einstein universe. They are used to investigate the 
>• ■ problem of how the curvature anisotropy affects the thermodynamic behaviors of an ideal 

^ ' gas. It is found that, when the universe is in the oblate configuration (i.e., d > 0), the 

effect of curvature anisotropy is to increase the number of the fraction in the Bose-Einstein 
condensation and to decrease the fermion distribution function at low temperature. When 
the universe is in the prolate configuration (i.e., d < 0), the effects of curvature anisotropy 
on the thermodynamic quantities is contrary to that in the oblate configuration. The den- 

O" 1 

sity matrix of a two particle system is evaluated and it is used to define the "statistical 
interparticle potential" as an attempt to give a "statistical interpretation" about the found 
' thermodynamic behaviors. It is found that when the universe is in the oblate (prolate) 

configuration the cnrvatnre anisotropy will enhance (rednce) botfi the "statistic attraction" 

S3- 



among the bosons and "statistical repulsion" among the fermions. It is expected that such 
a behavior will also be shown in the relativistic system. 
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1 Introduction 



The investigation of the thermodynamics of particle field in the curvature space-time is 
useful to understand the dynamics of the universe. Many efforts have been made to study 
the problem of Bose-Einstein condensation [1-3] in the static Einstein universe, [4-7] both 
in the nonrelativistic []4] and ultra-relativistic systems [5-7]. The condensation problem on 
other space-time such as the hyperbolic manifolds [8] or the Kaluza Klein universe [9] have 
also been studied. 

In the present article we will study some thermodynamic quantities of nonrelativistic 
ideal boson and fermion gases in the static Taub universe with a small anisotropy. Our 
main purpose is to see how the curvature anisotropy will affect the statistical properties 
of ideal boson and fermion gases. One basic reason for dealing with the static Tau space 
is that the complete eigenvalue spectrum of the Laplacian is known. This feature makes 
it mathematically tractable. There are several investigations of the physics in the static 
Taub universe [10], which include the quantum field effects [ll]and finite temperature field 
theory [12]. We will present the analytic forms of some interesting thermodynamic quantities 
expanded to the first order in a small anisotropy parameter d measuring the deformation 
from the spherical Einstein universe and thus several interesting physical properties can be 
read explicitly. 

The mode space-time and associated eigenfrequencies of a massless scalar particle [10] 
are briefly summarized in Sec. II. In Sec. Ill we study the phenomena of Bose-Einstein 
condensation along the algorithm of Al'taie [4] and find that, in the oblate configuration 
(i.e., d > 0) the effect of curvature anisotropy is to increase the fraction of condensation. 
In Sec. IV we show that the effect of curvature anisotropy (if d > 0) is to decrease the 
fermion distribution function at low temperature. The expressions derived in Sets. Ill and 
IV also show that effects of curvature anisotropy on the thermodynamic quantities of an 
ideal gas for the universe in the prolate configuration (i.e., d < 0) are contrary to that in 
the oblate configuration. As an attempt to give a "statistic interpretation" of those found 
thermodynamic behaviors we present in Sec. V an evaluation of the density matrix of a 
two-particle system and from this we define the "statistical interparticle potential" along 
the idea of Uhlenbeck [13]. We find that in the oblate (prolate) configuration the effect 
of curvature anisotropy is to enhance (reduce) both the "statistic attraction" among the 
bosons and "statistical repulsion" among fermions. Section VI is devoted to a discussion 
while some mathematic calculations are collected in the two Appendices. For simplicity 
we only study in this paper the case of a spinless particle, however, the case of spin-1 gas 
can be studied in a similar fashion. Note that although our investigation is confined to 
the nonrelativistic system, we expect that our main results of showing how the curvature 
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anisotropy effects the thermodynamic behaviors will also be found in the relativistic system 
[14-16]. For example, the Bose-Einstein condensation is a low-temperature phenomena and 
the system at low temperature can be described by the nonrelativistic effects [14]. However, 
to confirm this one can extend the analysis into the relativistic regime. This is more difficult 
in mathematics. We hope to be able to report on this elsewhere. 



2 Model 

The background geometry of the Taub space has the metric 

ds2 = -dt 2 + Y J V 2 M t )\ (2-1) 

i 

where i\ = £ 2 ^ £3 and the a 1 form a basis one-form on S* 3 satisfying the structure relations 
da' 1 = \<ji A a k . In the Euler-angle parametrization the a 1 are expressed as 

a 1 = cosipdO + sinil)sin6d(f). 

a 2 = —sintpdO + cosipsinOdcf). 

a 3 = dpsi + cos6d(f), (2.2) 

where < 9 < n,0 < ip, <fi < 2ir. The i\s are the three principle curvature radii of the 
homogeneous space and are constants for a static universe. The case when all three £\s are 
equal is the closed Friedmann- Robert son- Walker universe. The volume and curvature scalar 
of the Taub space is given by 

D _ 4ff-l§ _ 6(l + (4/3)rf 

U\ ~ a 2 (l + d) ' 1 j 

where a = 2£, and d = [£\j£%] — 1 is an anisotropy parameter with range — 1 < d < 00. We 
call the configuration with d > oblate and that with d < prolate. 

The eigenfrequencies of a massive (m) scalar field (4) with curvature coupling (£) to a 
static Taub universe are found to be [10] 

1/2 

(2.5) 




q 'U| ej K2 + m2 + l {1 - ()R 



Here L = J, K, M: J takes all values of positive integers and half integers and K,M = 
— J,—J + \,...,J — \,J. The quantum number M is totally degenerate. The nonrelativistic 
approximation of Eq. (2.5) can be expressed as a n-degeneracy frequency 

oo nq = -^[n 2 + d(n-l- 2q))% (2.6) 
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where n — 1,2, and q — 0,1, (n — 1). To obtain the above result we have discarded 
some overall constants which are irrelevant to our results. Notice that the eigenfrequency 
relation in the above equation, which is a nonrelativistic approximation and does not depend 
on the curvature coupling 6, can also be obtained by solving the nonrelativistic Schrodinger 
equation when the Laplacian operator is that on S 3 space defined by Eq. (2.2). Thus, we 
can regard Eq. (2.6) as the nonrelativistic approximation of the energy spectrum of an ideal 
gas, classical or quantum (boson or fermion), moving in the static Taub space and use it to 
investigate their thermodynamic behavior. 



Boson-Einstein Condensation 



The behavior of a nonrelativistic ideal boson gas confined to background geometry of an 
Einstein universe has been investigated by Al'taie [4]. Apart from the expected effects of 
smoothing out the singularities of the thermodynamic functions at the transition temperature 
he also verified a displacement of the specific heat maximum towards higher temperature and 
an enhancement of the number of condensed particles due to the finite size of the system. 
This section is presented to prove that such an effect of a finite size will be enhanced (reduced) 
by the curvature anisotropy for the Taub universe in the oblate (prolate) configuration. 
The number of particles for an ideal Boson gas in our specific case is 

oo n— 1 



-i 



n=l q=0 
oo n— 1 

= 12 n[exp{P'[n 2 + d(n - 1 - 2q) 2 ] + a} - l]' 1 

n=l q=0 
oo n— 1 

53 n[exp(p'n 2 + a) - = dn(n - 1 - 2q) 2 exp((3' n 2 + a){exp(/3'n 2 + a) - 1]~ 

n=l q=0 

d (_d_ d_ 
+ 3^ [dp' da, 



J2n 2 [expt((3'n 2 ) - 



(3.1) 



71=1 
1-1. 



where (3 = 1/kT and /j, = —(3 a is the chemical potential. We have expanded the result to 
the first order in the anisotropy parameter d and used the variables defined by 



(3' 



1 

47T 



A 

a 



A 



(—) 

\mkTj 



1/2 



2mkTa 2 

where A is the mean thermal wavelength of the particle. 

Substituting the Al'taie;s result [4] into the last summation in Eq. (3.1) we have 



(3.2) 



N = 



3^ \dp> da 



1/2 



r(-)F 3/2 (a)-27ra 1 / 2 S(2y) 



(3.3) 
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where 



= f>-< e— S(2y) = -^— i , y s 2^ g) . 



The function F^(a) has the properties (d/da)Fg(a) = —F^i(a) and 



'3/2 



c, (I) - 2(vra) 1 /2 + 0(a), F 1/2 (a) = it 1 ' 2 *' 1 ' 2 - t (\) + 0(a). 



(3.4) 



(3.5) 



At low temperatures a « 1, then using the above relations and after the calculation we 
obtain the result 



± c 



^-^ 2 + 7r2) (e^-l) 2 12. W 

(3.6) 

where £ is the mean interparticle distance. The term to the order of (A/a) 2 will be neglected, 
hereafter, because a >> A in our case. 

In the limit a — > oo and d = 0, Eq. (3.6) becomes 



(3.7) 



This is the well-known bulk result. Using the definition of the mean thermal wavelength of 
the particle in Eq. (3.2) we have the relation x = T/T = (X /X) 2 , then from the Eqs. (3.6) 
and (3.7) we get the final result 



1.3/2 

x) 



1 

1-- 

7T 



-2/3 £ ,-^1/2 

a \x. 



ycothy - - 



7T 2 e 2y 
—cothy - A(y 2 + tt 2 ) ^ _ ^ 2 ^ 



.(3.8) 



The case of d = is that in Ref. 4. As the anisotropy parameter d is small the "function " 
y can be expanded as 

y = y + dy 1 + 0(d 2 ). 
and Eq. (3.8) can be reduced into two equations 

~ 2/3 i(W /2 

a \xj 

.2 „2y 



y cothy = - 

7T 



[l-x-z' 2 ) 



(3.9) 



(3.10) 



., -/ 2 , 2^ e 



(3.11) 



Q(cothy - y csch 2 y ) 

For a given £/a and a given x Eq. (3.10) has been used to find yo by Al'taie [4]. Once the 
yo is found the value of y\ can easily be obtained from Eq. (3.11). Finally the value of y is 
then used to evaluate the thermodynamic quantities. 
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An interesting physical quantity is the condensate fraction N /N [1-3]. The number of 
particles in the ground in our case is given by 



N = [exp{(5' + a) - l]- 1 = 



r 

exp — 
W 



n -l 



+ 1 a 



4tt 3 



it 2 + y 2 a it 2 + y 2 y 2 



a 



4ir s x 
7i 2 + y 1 



Thus 



A^o _ / 1 A 2tix 
N ~ { + 2 ) n 2 + y 2 



5 J /a 



if 



3/2 



Using the relation (3.9) the above equation is reduced to 



N N y 7r 2 + ?/ 2 



7TX 



3/2 



|)rf[7T 2 + ^-% yi]. 



(3.12) 



(3.13) 



(3.14) 



To see how the curvature anisotropy will effect the Bose-Einstein condensation one must 
know the behavior of the bracket term in the above equation. Although we have as yet no 
ability to prove it, the numerical evaluation shows that the bracket term is positive definite. 
Thus we conclude that when the universe is in the oblate configuration (i.e., d > 0) the 
effect of curvature anisotropy is to increase the number of the fraction in the Bose-Einstein 
condensation. On the contrary, when the universe is in the prolate configuration (i.e., d < 0) 
the effect of curvature anisotropy is to decrease the number of the fraction in the Bose- 
Einstein condensation. This behavior may be interpreted as that the increase (decrease) in 
the "statistic attraction potential" among the boson gas by these effects, as shown in Sec. 
V, will push more (less) particles into the ground state. 

4 Fermion Distribution Function 

The total number of the fermion particles in our model is derived from 
N = Y / [exp[P(e-fi)] + l] 



oo n—1 

= ^^[ei#K-/ i ] + l]- 1 

n=l q=0 

oo n—1 

= EE n[exp{(3'[n 2 + d(n - 1 - 2qf] - £} + 1] 

n=l g=0 



-1 



oo n—1 

E n[exp([3'n 2 - f) + l]" 1 = dn(n - 1 - 2q) 2 exp([3'n 2 + a)[exp([3'n 2 - f) + 1]~ 2 

n=l q=0 
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d ol ( d <9 N 
1 + ^[w + dl / 



J2n 2 [exp(f3'n 2 -Z)-l}-\ 



(4.1) 



n=l 



where \x = f3~ 1 ^ is the chemical potential. We have expanded the result to the first order in 
the anisotropy parameter d and used the variables defined in Eq. (3.2). 

The last summation in Eq. (4.1) is executed in Appendix A and we have 



N 



1 + / lap** 



\ 3 /2 



1 7T 2 1 3 f(3 r 



1/2- 




1 7T 2 1 3 / d\ ((3 r 



1/2- 



(4.2) 



The above equation leads to the result 



N 
V 



6n 2 




(4.3) 



where g is a weight factor that arises from the internal structure of the fermion particle (such 
as the spin). The zeroth approximation of the above equation gives 



, TP 1 /6vr 2 iV\ 2/3 



(4.4) 



which is the ground-state result. Note that ej? is a function of particle density n = N/V 
which does not depend on the finite site radius, a, or on the curvature anisotropy, d. In the 
next approximation we obtain 



/i = kTl; = e F 



71 



'-»(£ 



1A fkT\ 1/2 
2 a \ne f ) 



1 + 



6. 



(4.5) 



which is the desired result. From the fermion distribution function n e = [1 + e^^- 1 ] -1 
and the above equation we see that the finite site effect (a < oo) and anisotropy effect 
(if d > 0) are to decrease the chemical potential and thus decrease the fermion distribution 
function. This behavior may be interpreted as follow: The increase in the "statistic repulsion 
potential" among the fermion gas by these effects, as shown in the next section, will spread 
particles into different states and thus decrease the fermion distribution function. 



5 Statistical Interparticle Potential 

It is well known that, in comparison with the normal statistical behavior, bosons exhibit 
a larger tendency of bunching together, i.e., a positive statistical correlation. In contrast, 
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fermions exhibit a negative statistical correlation. Uhlenbeck presented an interesting way of 
stating this property by introducing a "statistical interparticle potential" and then treating 
the particles classically [13]. This section will present an example to show how the curvature 
anisotropy effects the statistical interparticle potential, as an attempt to give a "statistical 
interpretation" of the results found in Sets. Ill and IV. 

Let us first briefly describe the method [3,13]. Define the one-particle matrix element of 
the Boltzmann factor by 

Fij =< Xile-^lXj >, (5.1) 

then the diagonal part of the matrix element of the Boltzmann factor for a system of two 
identical particles can be written as 



< Xi,X 2 \e H \Xi,X 2 >— FuF 22 ± Fi 2 F 2 i, 



(5.2) 



where the plus (minus) sign is adopted for the boson (fermion) system. For a translation 
symmetry system F 12 , is independent of the position and partition functions of a one-particle 
and two-particle system are 



Z l= Tr< Xi\e- pH \Xj >= VF tj = VF 

Z 2 =< X 1 ,X 2 \e- f>H \X 1 ,X 2 >= (VF ) 2 , 
and density matrix becomes 

< X ± , X 2 \p\X 1 , X 2 >= —{F\\F 22 ± F12-F21) = —r 

Z 2 V 



1 ± 



3. 

Fo 2 



(5.3) 



(5.4) 



The statistical potential v, is defined to be such that the Boltzmann factor exp(-flv) is 
precisely equal to the correlation factor (bracket term) in the above equation, i.e., 



v { r) = -kTin 



1 ± 



12 



(5.5) 



The problem reduces to show that the factor -^f- will be increased by the finite size factor 
and curvature anisotropy (if d > 0). 

Using the mode solution derived by Hu [10] we have 



00 n— 1 



F o = E E n exp{-(3'[n 2 + d(n - 1 - 2q) 2 ]} 



n=l q=0 
00 n—1 



E n exp(—/3'n 2 ) — nd{n — 1 — 2q) 2 (3'exp(—(3'n 2 ) 

71=1 q=0 
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1 + 5'( 1 + £ 



n 2 exp(—/3'n 2 ). 



(5.6) 



re=l 



The last summation in . (4.1) is executed in Appendix B and we have 

r 2 



dp 



7T 



1/4 



-(/?') 



/\-3/2 



(5.7) 



To evaluate F12 we consider the matrix element with the states ipi — 1^2 — 61 — O2 — <f>i — 
while 02 = 7T- (The other matrix element is difficult to calculate.) Using the mode solution 
derived by Hu [10] we have 



00 re— 1 

F 12 = X! X! n exp{—P'[n 2 + d(n — 1 — 2q) 2 ]}exp 
re=l <?=o 



(n - 1 - 2g)# 



00 j 

£ n exj9(-4^'K 2 )ea;p(-^ / n 2 )ea;p(i^) 

re=l JC=-J 



g ne -/3'« 2 [l_^'( n 2_ 1)]cos 



n=l 



2tt 



n-1) 



1 + + ^ 



n=l 



y(n-l) 



(5.8) 



The last summation in Eq. (4.1) is executed in Appendix B and we have 



F-\o — 



1 

= -e 



' - \ 3/2 1 



-1 



Thus 



(5.9) 



(5.10) 



which explicitly shows that the factor ^ will be increased by the finite size factor and 
curvature anisotropy (if d > 0). This means that in the oblate (prolate) configuration the 
effect of curvature anisotropy is to enhance (reduce) both the "statistic attraction" among 
the bosons and the "statistical repulsion" among fermions. 

We hope that further thermodynamic quantities, in addition to those derived in Sets. Ill 
and IV, can be interpreted by the statistical interparticle potential. 
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6 Discussions 



In the present article we study the phenomena of Bose-Einstein condensation in a static 
Taub universe and find that in the oblate configuration (i.e., d > 0) the effect of curvature 
anisotropy is to increase the fraction of condensation. We also present a possible way to 
see that the effect of curvature anisotropy (if d > 0) is to decrease the fermion distribu- 
tion function at low temperature. We have seen that effects of curvature anisotropy on the 
thermodynamic quantities of the ideal gas for the universe in the prolate configuration (i.e., 
d < 0) are contrary to that in the oblate configuration (i.e., d > 0). As an attempt to give 
a "statistic interpretation" of those found thermodynamic behaviors we have also evaluated 
the density matrix of a two-particle system and from this we defined the "statistical inter- 
particle potential" along the idea of Uhlenbeck [13]. We see that in the oblate (prolate) 
configuration the effect of curvature anisotropy is to enhance (reduce) both the "statistic 
attraction" among the bosons and "statistical repulsion " among the fermions. Although our 
investigation is confined to the nonrelativistic system, we expect that these properties will 
also be shown in the relativistic system. To confirm this property, a detailed analysis of the 
thermodynamics of a relativistic gas in the static Tau space-time is necessary. This is more 
difficult in mathematics and we hope to be able to report it elsewhere. 



APPENDIX A: Calculations of the summation in Eq.(4.1) 

Using the Poisson summation formula 



I OO />oo °° roc 

o/(0) + E /(») = / f(t)dt + 2j2 f(t)cos(2m7rt)dt, (Al) 



the sum in Eq. (4.1) can be written as 

i2 roc ^2 



71=1 L 



n 



oo j2 



+ e?' n2 -S Jo 1 + e^'* 2 -€ 



d,t 



m=l 1 



t 2 cos(2rmrt) 



d,t 



3/2 



I — 7 dx + 

Jo 2(l + e*-e + 



cos 



x 



1/2 



2mir I — J dx 



' 1 \ 3/2 
J') 



L + £ L r 



m=l 



(A2) 



For the case of £ >> 1, which is that at a low temperature, the integration in the above 
equation can be evaluated with the help of the generalized Sommerfeld's lemma [18] 



g{x) 



o 1 + e x -Z 



dx — I g(x)dx + 



T^fdi 
6 \dx 



+ 



360 I dx 3 



+ .... 



(A3) 
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Here g(x) is any well-behaved function. Thus 



7T 2 1 



(A4) 



m=l 
= g (/?') 3/2 



2m7r 



X 

J', 



1/2- 



m=l 



3tt3 



m J 7r 



m7r 



\ 1/2 



-j cos 



2™ (|) 



\ 1/2' 



2™ | 



1/21 



(A5) 



and we obtain the desired result 



3 U, 



3/2 



7T 2 1 3 



1 + 



1/21 



(A6) 



APPENDIX B: Calculations of the summation in Eqs.(5.6) and (5.8) 

Using the Poisson summation formula (Al) we have 



S = exp{—(3' n 2 ) 

n=l 

/•CO 00 ;-oo 

/ exp{-(3' t 2 )dt + 2 V / exp(-P' t 2 ) cos(2mirt)dt 
Jo m=l Jo 

1 Ar' 

KS 1/2(i+2e "^' 



i^-V /2 (i + 2 f>- m2 ^' 



m=l 



(5.1) 



which is a suitable approximation if a >> A [Note that j3' = (l/47r)(A/a).] The summation 
in Eq.(5.6) is calculated from 



n 2 exp(—f3' 



n 2 ) = 



OS 



9P' 

Using the Poisson summation formula (Al), the sum in Eq. (5.8) can be written as 



[B.2) 



n 2 exp{—(3' n 2 ) cos 



n=l 



7T 



n- 1) 
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y t e _/3 ' * 2 sin l^-tj dt + 2 J t e_/3 ' <2 sin [if) cos(2mnt)dt 



8 W 



3/2 



m=l 



7T 



+ (1 + 4m)ea;p 



1 /7T 

8 



W e -/3'- 2 /16/3' ( l_ 3e -W ); 



(S.3) 



which is a suitable approximation if a >> A. 
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